Conductance through a potential barrier embedded in a Luttinger Hquid: 
nonuniversal scahng at strong couphng. 
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We calculate the linear response conductance of electrons in a Luttinger liquid with arbitrary 
interaction g2, and subject to a potential barrier of arbitrary strength, as a function of temperature. 
We map the Hamiltonian in the basis of scattering states into an effective low energy Hamiltonian in 
current algebra form. First the renormalization group (RG) equation for weak interaction is derived 
in the current operator language both using the operator product expansion and the equation of 
motion method. To access the strong coupling regime, two methods of deducing the RG equation 
from perturbation theory, based on the scaling hypothesis and on the Callan-Symanzik formulation, 
are discussed. The important role of scale independent terms is emphasized. The latter depend on 
the regulaization scheme used (length versus temperature cutoff). Analyzing the perturbation theory 
in the fermionic representation, the diagrams contributing to the renormalization group /3-function 
are identified. A universal part of the /3-function is given by a ladder series and summed to all orders 
in g2- First non-universal corrections beyond the ladder series are discussed and are shown to differ 
from the exact solutions obtained within conformal field theory which use a different regularization 
scheme. The RG equation for the temperature dependent conductance is solved analytically. Our 
result agrees with known limiting cases. 



I. INTRODUCTION 

The quantum theory of charge transport in one- 
dimensional (ID) systems of interacting fermions is a fun- 
damental building block of our understanding of electron 
transport in nanostructures and of a future nanoelectron- 
ics. One-dimensional interacting fermion systems have 
been studied since the 1950s, beginning with the pio- 
neering works of Tomonaga and Luttinger. The early 
work exploited the fact that the elementary excitations in 
these systems are charge and spin excitations of bosonic 
character, and led to the formulation of the method of 
bosonizationjii^ This method offers a natural description 
of the fractionalization of the usual fermionic quasiparti- 
cles existing in higher dimensions, into spinon and holon 
quasiparticles in ID. While the bosonization method is 
highly successful in accounting for the properties of sys- 
tems in equilibrium, it is somewhat problematic when 
applied to transport situations. For example, an ideal 
quantum wire attached to two reservoirs is expected to 
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have a two-point conductance of Rq = ^ per channel. 
In contradiction to this general law it was found in early 
work on transport in clean Luttinger liquids^i^i^ that the 
conductance (here and in the following conductance will 
be measured in units of the conductance quantum Rq) is 
given by the Luttinger parameter K, and thus depends 
on the interaction strength. It was later pointed out that 
by carefully choosing the order of limits of frequency, 
— + 0, and wire length, L — > oo^^ or else by taking 
into account the screening of the external field by the in- 
teracting electron system^ one recovers the correct value 
of unity for the two-point conductance. Whether the for- 
mer or the latter explanation is the correct one remains 
disputed to this day. 



The more relevant case of a Luttinger liquid with po- 
tential barrier has been considered first by Kane and 
Fisher=^ and by Furusaki and Nagaosa,— again using the 
bosonization method. These authors found that interac- 
tion has a dramatic effect on the conductance: for repul- 
sive interaction the conductance is found to tend to zero 
as a fractional power of the temperature T , in the limit 
T — !■ . This was shown in the two limits of a weakly 
scattering barrier and a strong (tunneling) barrier. In 
addition, results in the complete temperature range are 
available at special values of the interaction parameter, 
K = \ and K = 1 ^i^iiaiiiii^ All these works suffer from 
the drawback that in the clean limit the conductance 
is found to tend to K rather than 1 . It has been ar- 
gued that these results may be applied to the four-point 
conductance. However, the four-point conductance is ex- 
pected to tend to infinity in the limit of vanishing barrier 
strength, a behavior not shown. 

For the reasons noted above we think it worthwhile 
to develop an alternative formulation of the transport 
theory of Luttinger liquids, formulated entirely in the 
fermionic language. The fermionic representation offers 
the possibility to connect the fermionic degrees of free- 
dom in the (non-interacting) leads smoothly with the 
degrees of freedom in the interacting system. In other 
words, it allows to use the scattering states of the sys- 
tem in the non-interacting limit as a basis of description. 
On the simplest level, in lowest order in the interaction, 
this program has been carried out in a seminal paper by 
Yue, Matveev and Glazmanii^ The physics of this prob- 
lem lies in the scale-dependent build-up of a polarization 
potential around the bare barrier, induced by the Friedel 
oscillations of the density. For repulsive interaction, the 
polarization potential is found to extend farther and far- 
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ther out as the temperature is lowered, until at T = 
it is infinitely extended, leading to a vanishing transmis- 
sion probability across the barrier. The process of grad- 
ual growth of the effective potential barrier may be de- 
scribed in the language of renormalization group theory 
applied to the transmission probability as a function of 
the temperature. A generalization of the approach of Yue 
et al. to the case of two barriers has been given in Ref. 
p^ . These predictions, based on the continuum version 
of the theory, were thoroughly checked and confirmed by 
the fermionic functional renormalization group method, 
starting from the Hubbard models on a latticeJ^ii^ 

In this paper we report on a significant extension of 
the work of Yue et al. to the case of systems of spin- 
less fermions in Id interacting via arbitrarily strong for- 
ward scattering (parameter (72) Siiid subject to a short 
range barrier potential (width a). The principal tool we 
will be using is perturbation theory in the interaction, 
summed to infinite order. To achieve this goal, and to 
gain insight into the structure of perturbation theory, 
we map the problem first onto an effective model of in- 
teracting currents of chiral fermions. The reformulation 
allows to describe the effect of the barrier potential as 
a local magnetic field at the position of the barrier, act- 
ing on the pscudospin vector of the chiral current. The 
logarithmic corrections to this field, stemming from the 
fermionic interactions can be addressed in a simplified 
poor-man scaling approach; however this approach be- 
comes ambiguous in higher orders. Therefore at a later 
stage we will apply the Callan-Szymanzik type of anal- 
ysis, by first introducing the counter terms to magnetic 
field part of the Hamiltonian to compensate the effect of 
the interaction, and secondly requiring the "bare" value 
of the field to be independent of the high-energy cut- 
off. Equivalently, we may assume the scaling property 
of the conductance to hold, entailing the existence of a 
renormalization group equation for the conductance as a 
function of the scaling variable A (at finite temperature 
A = ln(To/r) ). By comparison with the structure of per- 
turbation theory for the conductance, which is given by a 
power series in the interaction and in the scaling param- 
eter A , one may extract the RG /3-function. In this way 
we derive a universal renormalization group equation for 
the conductance as a function of the scaling parameter. 

The most important consequence of the current alge- 
bra formulation is, however, to allow for an efficient or- 
ganization of perturbation theory. After presenting the 
formulation of perturbation theory as well as a few low 
order results we analyze the structure of the perturba- 
tion expansion with respect to logarithmically divergent 
terms containing powers of A = ln(_L/a) or ln{To/T). 
We identify a class of universal diagrams contributing the 
principal terms linear in A and derive an integral equa- 
tion resumming these contributions ( "ladder approxima- 
tion"). These terms are shown to dominate in both lim- 
its, small and large conductance. We identify the pref- 
actor of the terms in the perturbation series for the con- 
ductance G linear in A with the renormalization group 



/3-function of the flow equation of G(A). There is, how- 
ever, one important new feature: there appear scale in- 
dependent terms (from third order on), which lead to a 
correction of the conductance even at A = 0, the ultra- 
violet cutoff. These terms must be taken into account, 
in order to preserve the scaling property of the conduc- 
tance, and lead to a starting value of the conductance 
G{A = 0), different from the conductance in the absence 
of interaction. 

The RG equation may be solved analytically in the lad- 
der approximation, to give the conductance as a function 
of temperature for any interaction strength and barrier 
potential. The result completely agrees with earlier work, 
where a comparison is appropriate We emphasize that 
our result for the conductance reduces to unity in the 
absence of a barrier. Furthermore, the effect of screening 
of the external field'^ is included by omitting all polar- 
ization diagrams. We checked the renormalizability of 
the theory explicitly by calculating all terms up to third 
order in A and 52 ■ 

At intermediate temperatures (semi-transparent bar- 
rier) additional diagrams contribute small corrections to 
the /3-function. We calculate these terms in lowest (third) 
order in (72 ■ A proper treatment of these contributions 
requires taking into account scale independent terms in a 
way sketched above. We demonstrate explicitly that the 
resulting /3-function is uniquely determined for a given 
physical quantity. However, (3 turns out to be a slightly 
different function for the two cases of interest here, the 
temperature dependent and the length dependent con- 
ductance. The result of the solution of the RG-equation 
for the conductance, using the approximate /9-function 
thus obtained, is found to agree quantitatively with all 
known results on the scaling behavior of the conductance 
in the limits G — *■ and G —^ 1 . In the intermediate 
regime G ^ 1/2 we find small discrepancies with the re- 
sults of conformal field theory in the case K = ^, which 
we can trace back to a different cutoff scheme used there 
and the effect of higher order terms neglected in our work. 
As an excellent interpolation formula between two exact 
limits our result goes beyond all previous works: it is 
valid for any interaction strength K and any potential 
scattering strength. 



A. Formulation of the problem 

The Hamiltonian includes the kinetic energy of freely 
right- and left-moving (R,L) spinless fermions with lin- 
earized dispersion, Hq, the interaction energy between 
the left- and right-moving densities, iJi, and the impu- 
rity part decribing scattering off the barrier, placed at 
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the origin, Himp : 

H = Hq + Hi + Hi, 



Hamiltonian Himp in ([1} is then written as 



(1) 



Ho = vf 



dx 



Hi = .92 / dx {■>PI^'iPb.){'>PI'4'l) 



In order to regularize the theory, we assume that the 
region of interaction |a;| > a, where 52 7^ 0, and the region 
over which the impurity potential is nonzero, |a;| < a ^ 
kp (kpis the Fermi wave vector) do not overlap spatially. 
The length scale a will serve as an ultraviolet cutoff in 
the scale dependent logarithmic corrections considered 
below. 

The last term, i?impj in GI is not explicitly specified 
here (but see below). We will rather take the one-particle 
scattering states to be known. This means that we char- 
acterize the barrier by transmission and reflection ampli- 
tudes t = i = cos 9, r = — f* = zsin^e"'' , with negligible 
energy dependence in the energy range of interest (here 
t, t, r, f are the usual components of the impurity scat- 
tering matrix). Notice also that we assume the system 
to be non-interacting in the leads, > L, which allows 
for an asymptotic single-particle scattering states repre- 
sentation. 

We define the creation operator of scattering states as 



^C2fe, 2; < 



V'i.(x) = (e^'=--^rfce-^'=-)4-|-4e-^'=-4 

= tfee^'^^cl^ -I- (ffee^'=" -t- e-^^-")4fe, x > (2) 



with c\f, (c\f^) the creation operators of asymptotically 
right-going (left-going) fermions of momentum k. Such 
a representation requires the kinetic energy part of the 
Hamiltonian to be of the form — V^/(2m), and the mo- 
mentum, k = \/1mE > 0, before the linearization of 
the dispersion around the two Fermi points. In appendix 
A we clarify the correspondence between the scattering 
states representation and our linearized Hamiltonian 
([1]) written in the basis of chiral fermions. 

For simplicity we do not consider the so-called part 
of the fermionic interaction in this work, i.e. the term 

j'^^dx [{tplipRf + ii^li^L^]. For L ^ cx) the 54- 
interaction can be absorbed into the redefinition of the 
Fermi velocity, vp = vp + g4/2Tr. For finite L one can 
show that g4 does not lead to a renormalization of d.c. 
conductance, which is the quantity of interest below. The 
effect of g^ on the a.c. conductance can be analyzed, 
following the guidelines im^^. 

It appears that in the case of weak potential scatter- 
ing, U{x) <C Ep, and for electrons with energies close to 
the Fermi energy, Ep, the situation is characterized two 
limits of the Fourier transform U (g) of the potential: the 
forward scattering amplitude, U{q ~ 0) = Ui, and the 
backward scattering amplitude, U{q 2kp) = U2, with 
real- valued ui and complex- valued 1*2- The fermionic 



H, 



imp 



— Vp I dx 



ui(x)(^JjVfl, + ^l^L) (3) 



+ {u2{x)4'\i^R. + h.C. 



where the dimensionless functions ^1,2(2;) are short-range 
impurity potentials, which means that ^1,2(2;) — for 
\x\ > a and the above amplitudes ui^2 = J dxui^2{x). 

In Appendix Rl we show that the connection between 
the microscopic Hamiltonian and the 5'-matrix can be 
clarified in some simple cases, but requires a detailed 
knowledge of ui^2{x), when one goes beyond the lowest 
Born approximation. 

Returning to we define Fourier transforms, 

ikx ~-h 



V-^" (^) = r He^'^St. and i^tixl = f e- 
Then we have for the electron creation operator at posi- 
tion X 

ij+{x) = {e{-x)[2P+{x)+ri;+{-x)+i^4{x)] 

+ Six) [t^+{x) + Pi^+i-x) + i;+{x)] } , (4) 

with the step function 9(a;) = 1 at x > 0. 



II. INTERACTING CASE, PREVIOUS WORKS 
A. Boundary sine-Gordon model 

In bosonization technique, the chiral fermions are rep- 
resented as exponentials of chiral boson fields, ipR ^ 
f,-iVR^ ~ e*'^^. with the fields ipR(^L) ^ V T 0. 
Here the primary field and its canonically conjugate 
momentum H — ir^^dxO obey the commutation rela- 
tion [(^(a;), H(j/)] — iS{x — y). We have for the density 
'4'r'4'r + i'L^'L — dxip/TT and the Hamiltonian ^ 
can be represented as 



H 



2tt 



dx 



K {dxOf + K'\dxvf 



+2ui{x)dLp{x) + U2{x) cos2ip{x) 



(5) 



One has vp = vp and K ^ I for free fermions, whereas 
the short-range interaction Hi between the fermions, Eq. 
(HI), leads to the renormalized Fermi velocity vp = vp{l — 
g^Y^"^ and the Luttinger parameter K = [(1 — 5)/(l -I- 
(7)]^/^. Here and below we use 



a = g2l{'2i^vp) 



(6) 



One usually argues, that the term ui can be absorbed 
into a redefinition of ^^(x) — > ^^(x) -\- uisgn{x) (cf. Ap- 
pendix After this redefinition one concentrates on 
the U2 term which constitutes the boundary sine-Gordon 
(BSG) problem. A solution to this problem, in the lowest 
order of U2 was presented in the early work.^'^ Focussing 
on the small-energy sector of the problem, one uses the 
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renormalization group approach, removing the higher en- 
ergy states of the problem while simultaneously rescaling 
the parameters of the action. It turns out that under this 
procedure the amplitude of the BSG term is renormalized 
as 

U2 ^ U2e(i-^)^ (7) 

where the RG scaling variable A = hiEp/e > . As a 
result, the conductance at small U2 , given by G = |tp ~ 
1 — , undergoes a similar renormalization 



G~ 1 



I |2 2(l-if)A 



(8) 



with the lower cutoff in A defined by temperature or volt- 
age bias, A = ln£^;-/ max[T, V^]. We shall confine our- 
selves to the the linear response regime, V ^ T in the 
following. 

For repulsive interaction K < 1 the conductance de- 
creases with falling temperature. When the renormalized 
impurity potential becomes strong, |M2pe^^^~^^^ ^ 1, 
the weak-impurity assumption is violated and one should 
use a different line of reasoning. One may think of two 
semi-infinite parts of the wire, connected by a weak tun- 
neling link t. An appropriate RG treatment shows that 
the repulsive interaction leads to a further decrease of the 
conductance, which acquires the form 



G 



^l)A 



(9) 

These two regimes, Eqs. ([5]), © show different scal- 
ing exponents, which are approximately equal for small 
interaction 



1- K ~K- 



A smooth crossover between the two asymptotes ([5]) and 
^ was first explicitly provided for a special case of in- 
teraction strength, K = 1/2, which is known as Luther- 
Emery refermionization point. The full curve for the 
conductance was obtained^iiP in the form of a one- 
parameter scaling function. 

In subsequent works^^'^^ the BSG model was analyzed 
at arbitrary values of Luttinger parameter, K. The set 
of coupled integral equations was shown to be finite for 
1/K = 1,2,3. . . and the solution for = 1/3 was de- 
rived, in particular. Corresponding scaling functions for 
the conductance were obtained for finite T and voltage 
bias on the contact V either by numerical solution of the 
corresponding integral equation or as a series in powers of 
the scaling parameter, see also^. For practical purposes 
this form of representation is not very convenient. 

In principle, the statement of a one-parameter scal- 
ing function for the conductance and the knowledge 
of its actual form provides a solution for the temper- 
ature renormalization of the barrier apparently at ar- 
bitrary strength. However, it was noted iniiii^, that 
"the strong-barrier problem which is at the end of our 
renormalization-group trajectory follows formally from 
dimensional continuation of the integrals for the weak 
barrier problem. It is not in any case a generic strong- 
barrier problem." 



B. Fermionic approach, RG for S-matrix 

Yue, Matveev and GlazmanAS have developed a theory, 
starting from the formalism of scattering states and tak- 
ing into account the fermionic interaction in lowest order 
of perturbation. They arrived at the RG equation for the 
transmission amplitude t in the form 



dt , , ,,, 

^ = -.t(l-|a 



The solution to this equation is found as 



t = 



to 



(|toP + (l-|ioP)e2^"^)i/2 



(10) 



(11) 



which, in leading order in g, agrees with Eqs. ([H), ([9]). 

The advantage of this approach is its applicability for 
arbitrary strength of impurity potential. One may thus 
avoid the step involving the transition from the initial 
Hamiltonian ^ to the observed conductance, which as 
we saw above can depend on the short-distance (ultravi- 
olet) details of relevant quantities. This view is further 
corroborated by the work by Callan et ali^ who dis- 
cussed the solution of the non-interacting {K = 1) BSG 
theory, and arrived at a current algebra formulation sim- 
ilar to our formulation below. These authors show that 
the connection of the initial model ^ to the observable 
quantities depends on the ultraviolet regularization. We 
will return to this point in our discussion below. 



III. SCATTERING STATES AND CURRENT 
ALGEBRA 

In the above we defined a scattering states representa- 
tion (HI in chiral representation. In order to describe the 
effects of interaction, we need to define fermionic densi- 
ties in the same basis. One may form four bilinear density 
combinations out of two chiral fermions. It is convenient 
to group these combinations into chiral densities or "cur- 
rents" , according to 



J{x) 



tpl{-x)lpi{x) ljjl(-x)lp2{-x) 



Jo 

Jl- 



iJ2 



Ji - iJ2 
Jo — J3 



(12) 
(13) 



We call Jo the pseudocharge current and the vector J ~ 
(Ji, J2, J3) the pseudospin current. These operators obey 
U (1) and SU (2) Kac-Moody algebras, respectively^^, as 
we discuss below. 

The particle density operators for incoming (i) and 
outgoing (o) particles in terms of the J^'s are given by 
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(here and in the following x > ) 

P-mi-x) = i1jI{~x)iJji{-x) ^ Jo{-x) + Jsi-x) (14) 

PiL{x) = ■il;2{x)'ip2{x) ^ Jo{-x) - J3{-x) (15) 

PoBix) = {ti;\{x)+f^jl{-x)){t*Mx)+f-*M-x)) 

= {S.J{x).S^)n^Mx) + J3{x) (16) 

Pol{-x) = (rV4(x)+tV4(_2.))(^*^^(3,)+f*y,2(„3,)) 

= {S.J{x).S^h2^Jo{x)-Mx) (17) 

Here J3 = {RJ)3 is the third component of the pseu- 
dospin current vector J rotated by the orthogonal matrix 
Rfiu given by 



Riii^ = ^Tr[(Tf,.S.a,,.S'^] 



(18) 



with Pauli matrices cr^. It is seen from ([T51) that the 
global phase of S drops out completely. Omitting this 
phase, we parametrize 



S ■ 



t r 
r t 



cos f 



I smfe 
cos^ 



(19) 



In this notation, the components of interest below read 
i?33 = W - = cos 261, i?32 = Im{<P + tr*} = 
— sin 26* cos 0, and R31 — Re{tr* — tr*} = sin 26* sine/). 

We see that the description of the barrier is given by 
the eigenmodes of the current Ji{x), with negative and 
positive x and i = 0, . . . , 3. In this picture, the incom- 
ing current is given by the diagonal components of the 
J matrix, and the outgoing currents correspond to the 
diagonal components of the rotated matrix S.J.S^. Ev- 
idently, the pseudocharge component Jo is not affected 
by this rotation. 



A. Formal identities and Green's functions 

The chiral densities introduced above satisfy the Kac- 
Moody commutation rules 



[Jo{x),Jo{y)] = —dx6{x-y) 



(20) 



[Jj{x),Jkiy)] = -^SjkdxS{x-y) + iSjkiJi{y)S{x-y) 

with j,k,l — 1, ... 3 and totally antisymmetric ejki- The 
short-distance operator product expansions (OPEs) are 



Mx)Jo{y) 



-1 1 



Svr^ (.T - y - iOy 



T f \ T f \ , ^ £]kiJi{y) 

Svr-^ (x — y — lO)^ 2n x — y — lU 

In simple terms, relevant to our discussion below, one can 
consider (|^D|) as a consequence of (HD). In their turn, the 
OPEs (j2ip simply represent fermionic Green's functions, 



[27r(a; — y)]~^- A bilinear form in the currents is a prod- 
uct of four fermion operators; the complete convolution 
of these gives a product of two two-point Green's func- 
tion (modulo normally ordered operators), and the par- 
tial convolution produces a Green's function multiplied 
by bilinear forms in the fermion operators, i.e. current 
operators. The position argument of the latter current 
operators (last term in (HI])) is subject to ambiguity, but 
usually this concerns only less relevant terms. In our 
analysis below we deal with an action non-local in the 
currents, and to handle such an ambiguity becomes an 
increasingly difficult matter. This is why we resort to the 
perturbative treatment in terms of chiral fermions. 



B. Observable current and conductance 

The physical charge density is p{x) = PiR(x) + Pol{x) 
at a; < and p{x) = Por{x) -\- PiL(x) at a; > 0. A simple 
calculation shows that 

P{x) = Pc{x) + Ps{x), 

Pc{x) = Jo(-a;)-f Jo(a;), (22) 
Ps{x) = sign{x) {-,h{-\x\) + J3{\x\)^ . 

The physical current operator is also decomposed into 
two parts, in the pseudocharge and pseudospin sectors. 
We use the continuity equation —dxj{x) = dtp{x) = 
—i[p{x),H] to obtain 



i{x) = ic{x) +js{x), 
jc{x) = vf{-Jo{-x) + Jq{x)), 
3s{x) = vf{J3{-\x\) + J3{\x\)), 



(23) 



It follows then that the pseudocharge current corre- 
sponds to the part of physical density which is even with 
respect to the reflection at the scatterer, pc{—x) = Pc{x). 
The pseudospin current corresponds to the odd part of 
the density, ps{—x) = —ps{x). This symmetry prop- 
erty simplifies the discussion of a voltage bias, applied 
to the barrier. One sees immediately, that applying the 
same electrochemical potential to both leads of the wire 
corresponds to its coupling to the pseudocharge current. 
It results in an overall increase of the fermionic density, 
symmetric with respect to the barrier, and not to a net 
physical current. By contrast, an electric potential of the 
form V{x) = }-V sign{x), couples only to the pseudospin 
part, leading to a term in the Hamiltonian 

Hv = V dx {J3{x) - M-x)) 

JQ 

The physical current induced by a bias voltage is there- 
fore given in linear response theory as 

{j{x,Lu)) = Gix,uj)V{Lo) 

G{x,t) = ^iQ{t){[js{x,t), / dyps{y,0)]) (24) 

Jo 
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FIG. 1: (Color online) (a)The interaction of the left- and 
right-going densities in the basis of scattered waves, (b) The 
interaction of the left- and right-going densities in the chiral 
basis, corresponding to non-local interaction, Eq. (|26|) . 

A remark is in order here. Strictly speaking, the prefactor 
vp in the definition (|23p applies for non-interacting leads. 
In the interacting region we have to replace vp vp — 
(72/27r in (|23p . We may avoid this complication, analyzed 
elsewhereiiii, by assuming in ([M]) that both js{x, t) and 
Ps{y,0) are placed outside the interacting region, x,y > 
L. The change of the limits of integration over y in (124]) 
is not important in the limit of d.c. conductance \uj\ <C 
vp/L, considered below. 

C. Hamiltonian 

It is known that fermions with linearized dispersion in 
one spatial dimension can be described by a Hamilto- 
nian, quadratic in chiral fermionic densities. This obser- 
vation, usually attributed to Tomonaga, is the basis of 
the bosonization approach.^ Similarly to ([5]), we have 

/>oo 

Jo 

which can be represented as 

Ho = 2ttvf / dx {Jl{x) + J'i^x)) (25) 

J —OO 

poo 

+27rvp / dx {Jq{x) + Jl{x)) 
Jo 

The interaction part is 

Hi = §2 i dx{piR{-x)poLi-x) + PoR{x)piL{x)) 

J a 

= 2g2 [ dxiJoi-x)Joix)~M-x)J3{x)) (26) 

J a 

In Fig. [T] we show in a pictorial way our parametrization 
of the fermionic densities and the interaction. Hi. In Fig. 



[T^ the (72-iiiteraction processes are shown in the usual 
scattering configuration. The representation in terms of 
the chiral currents (all particles moving to the right) 
is shown in Fig. [TJa and leads to a seemingly nonlocal 
interaction. Notice that in the case of perfect reflection, 
t = 0, we have J3 — — J3 and the observable densities 
(see below) form an Abelian U{1) sub-algebra of SU{2). 
This case of "open boundary bosonization" , considered 
by Fabrizio and Gogolin, allows a complete and rather 
simple analysis^i . In the next subsection we elucidate 
what happens at the origin, x — 0. . 

D. The barrier as magnetic field 

Using the parametrization (|19p. one finds that Rij = 
exp(£yfci?fe) i.e., it is a finite rotation characterized by 
the dual vector B = 26'(cos 0, sin 0). We cannot sym- 
metrize the Hamiltonian in the pseudospin space, reduc- 
ing it to the so-called Sugawara form,— because the vector 
B breaks the SU (2) symmetry of the problem. We may 
symmetrize the Hamiltonian in the plane perpendicular 
to B, but that provides little advantage . 

The above formulation of the problem, incorporating 
the knowledge about the scatterer in the asymptotic scat- 
tering states, is equivalent to the the alternative formu- 
lation of the pseudospin current being rotated by a ficti- 
tious magnetic field B located at the origin. The pseu- 
docharge sector remains unchanged and the pseudospin 
sector of the problem can be equivalently described by 
the following Hamiltonian 

/OO 
dx {2nvFOl{x) - g203{-x)03{x) 
-00 

+vfBjOj{x)S{x)) , (27) 

where the operators Oj{x) = U^Jj{x)U, with j = 1,2,3 
obey Kac-Moody relations (PD|) . The two Hamiltonians, 
(PS)) . and (PT)) . are connected by a canonical trans- 
formation H' = U^HU , where 

/•OO 

U — exp i / dxBj Jj (x) . 
Jo 

We derive the relation between Oj{x) and Jj{x) in Ap- 
pendix |B] 

Equation ([7f|) resembles the Hamiltonian for the 
Kondo problem in the current algebra representation!^ 
However there are important differences. The first one 
is the non-local interaction (72 between the incoming and 
outgoing waves in our case. The second difference is the 
classical nature of the field B, as opposed to the quantum 
nature of the Kondo spin S. 

The matrix Green's function for chiral fermions is di- 
agonal in pseudospin space when defined in terms of the 
asymptotic eigenstates. However in terms of the usual 
right- and left-movers its matrix structure is evidently 
given by Eq. (fTO|) . which shows the action of the mag- 
netic field S = expisB, with s = ^<t. We will use this 
fact in Sec. I VII below. 
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IV. RG EQUATION FOR 5-MATRIX 

Let us show how the renormahzation group equation 
for the many-body S*— matrix is simply reproduced in 
the current operator formahsm. FoUowing Affleck and 
Ludwig,^" we write the 52 term of the 5*— matrix in the 
interaction representation as 

Ttexp2ig2 J '^^ J dxJ3{—x){RJ)3{x) 

Expanding the latter object in first power of 52, we find 
at given t 



dxJ3{—x)R3jJj{x) 



Applying OPE rules (|2T|) to the last expression one finds 
an operator term 



~2ti'^2x 



After eliminating the distances close to the scatterer, i.e. 
integrating over x from a to a certain a', which serves 
now as new lower cutoff, we obtain the term which can 
be re-exponentiated to the action giving a contribution 
of the form SH = vf SBk Jfc(O), with 



3j 



92 



dx 

X 



We notice further that e^jkRsj = Bk — 
sin 26* (cos sin 0, 0), i.e., the additional magnetic 
field at the origin is parallel to the above B, and we 
may restrict the discussion to the change in modulus, 
B ^ B + 6B. Denoting dA = J" dx/x, we obtain 
d(26')/dA = gsin26i, or 



d sin^ e 
dA 



= 2g siiT e cos^ 9 



(28) 



which is equivalent to Eq. pO)) . 

After this verification of the previous result in the first 
order in we are tempted to obtain the next-order cor- 
rections to the RG equation It turns out however 
that working entirely in the current operator formalism 
gives no advantage in such analysis, as we explain in the 
next subsection. 



A. Higher orders in current operator formalism 

Let us return to the last term in (PT|) 



Using the definition (fT2l) . and free fermionic Green's func- 
tions, one obtains e.g. the precise relation 



J2{x)J:i{v) 



87r(a; - y) 



(29) 



Jo{x)Jk{y) 



2ti{x - y)' 



modulo the normal ordered part. Evidently, the numer- 
ator in (j29p is reduced to 4Ji(y) only in the limit x = y. 
In other words, the relation ([2T|) is the "short-distance" 
expansion and must be used with caution. The meaning 
of the relation l|2ip becomes clearer when we compare 
the results obtained within the purely fermionic (and 
thus unambiguous) approach with the results provided 
by the use of (HI]). When available, this comparison 
yields identical results in the limit of large (external) 
times and distances. Hence, the relations (I^Tj) provide a 
good mnemonic rule for a quick estimate of a given con- 
tribution to a desired quantity. However, if we want to 
obtain concrete prefactors and verify the renormalizabil- 
ity of the theory, it is better to use the standard fermionic 
approach, treating the currents Ji as pseudospin vertices. 
Proceeding this way below, we use a symmetry of the 
problem, namely the unitarity of the S'-matrix related to 
the conservation of charge, which leads to a pure rotation 
of the vector current operator after the scattering event. 

A different attempt to explore the current structure 
of the theory would be the use of the equation of mo- 
tion method, which we describe in the Appendices [Cl and 
Id] If successful, this strategy would account for the in- 
fluence of interaction in all orders. The main difflculty 
arising in this approach is the nonlocal character of in- 
teraction (|26p and the appearance of an infinite series of 
coupled equations. We show in Appendix [D] that the 
attempt to truncate this series by considering the " most 
relevant" contributions leads, perhaps not surprisingly, 
to the restoration of the previous result (fTT]l . In order to 
go beyond this level of consideration, we use the regular 
perturbation theory in chiral fermions in the remainder 
of the paper. 



V. CORRECTIONS TO CONDUCTANCE: 
PERTURBATION THEORY 

A. Diagrammatic technique 

Let us first formulate the rules for the calculation of 
corrections to the conductance in perturbation theory in 
the interaction g. 

The contributions to G{flm) in n-th order of g2 may 
be calculated with the help of Feynman diagrams in the 
position-energy representation (ri,„ is the external Mat- 
subara frequency). We draw n vertical wavy lines in par- 
allel, each carrying the factor —2g2, the upper endpoint 
of the i-th line at —Xi with pseudospin matrix ^t^^ , the 

lower one at Xi with matrix 5^3/^1"^^ attached ; a, {(3) are 
pseudospin indices of incoming (outgoing) fermion lines. 
The external vertices are at —y with matrix and at 
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X with matrix ^{Rf)^. The vertex points are connected 
by Green's functions 

g{x ; LJ„) = -ii'^'signK)eKa;)e-'"""/''^ , (30) 

where the uj„ are Matsubara fequencies uj„ = (2ri + l)7rT. 
All internal x -variables are integrated on the positive 
semi axis. The trace over the product of all isospin ma- 
trices in each fermion loop is taken and a factor of l/nl is 
applied to each n-th order diagram. The limit fl„i +0 
is taken at the end. 

After this brief overview we provide further details 
about our calculation. 

(i) From Eqs. ((221), (El) it follows that only two out 
of four correlation functions in G{x,uj) are non zero. 
Namely, the chiral character of the currents (PT|) leads 
to {J3{x,t)J3{y,0))uj = at X < y. Then assuming 
y > X = L in ([24|) we get 

G(x,t)~ J dy{[M-x,t) + Mx,t),M-y,0)]) 

The first term here gives a contribution G''-^-'(a; ^ 0) = 
1/2, it is not affected by interaction. The second term is 
G(2)(w ^ 0) = R33/2 = cos(26')/2 for free electrons, so 
that G = G(i) + G(2) = cos2 9 = jtp as expected. 

(ii) It is this second correlation function, ~ 
( J3(a;, t) J3(— y, 0)), which undergoes renormalization. 
Let us denote this quantity by Y : 

y = 2G - 1 = i?33 = cos 261 (31) 

The decrease of the conductance can be viewed as the 
increase of the rotation angle 9. Below we discuss the 
corrections to Y. 

(iii) In order 5", we draw a diagram which contains 
2n -|- 2 current vertices. Two currents are external, and 
2n come from n pairs of currents in each interaction term. 
We fix the coordinates of the interaction terms and cal- 
culate the diagrams. In the limit of external frequency 
r2 — > each diagram is cx 51™ with m > 1, the integration 
over external y in (I24p reduces this power m —^ m — 1 
so that only the linear-in-r2 contribution survives in this 
limit ; see also paragraph (viii) below. After this pro- 
cedure we have finite contributions, dependent on the 
positions of interaction points, Xi, and should integrate 
over these positions. 

(iv) These last integrals over Xi (i = 1, . . . n), if taken 
over the entire semiaxis (0,cxd), diverge logarithmically 
and it is precisely at this step when we have to introduce 
the cutoff, as we discuss at length below. 

(v) Each current vertex has matrix structure and is 
(T3 for a current at negative coordinate (before scatter- 
ing) and cos 29as + sin 29a2 for positive coordinate (after 
scattering). We denote these matrices by Vj. Correc- 
tions are generated by connecting all current vertices by 
Green's functions. Technically, we employ Mathematica 
and generate all possible permutations {ji, j2, ■ • ■ ,j2n+2} 
of vertices and make one fermionic loop, connecting these 



vertices in given order. The advantage of current repre- 
sentation is that the information about the scattering is 
encoded in Pauli matrices and the kinetics is given by 
chiral Green's functions pop . The matrix structure pro- 
duces an overall prefactor before the diagram of the form 

lTr[Vn-Vn---- 

(vi) Other corrections are generated when 2n + 2 cur- 
rent vertices are connected by more than one fermionic 
loop. We denote these structurally different sets of dia- 
grams by listing the number of current vertices in internal 
loops. In the zeroth order of g we have the trivial set {2}, 
in the first order of g we have only a set {4}. In second 
order of g we have two sets, {6} and {4, 2}. In the third 
order of g we have four different sets : {8}, {6, 2}, {4, 4} 
and {4,2,2}. Some of these diagrams are shown in the 
Figures [21 and El 

(vii) To exclude double counting of different graphs in 
this procedure, we should fix one element in each loop. 
We always fix one external vertex as a first element in the 
first loop. In addition, we assume that the coordinates of 
internal vertices are ordered, i.e. L > xi > X2 > X3 > 0, 
etc. Then. e.g. the number of graphs in the {8} set be- 
comes 7! = 5040. In the {6,2} set, one has 7!/2 graphs, 
with factor 2 stemming from the necessity to fix one el- 
ement in the second loop. Similarly, the set {4,4} is ob- 
tained by taking 7! permutations, then partitioning each 
permutation into two parts of length 4, summing all pos- 
sible diagrammatic contributions and dividing the result 
by factor 4 to account for the necessity to fix one element 
in the second loop. 

(viii) We do not explicitly require that both external 
vertices belong to the same loop, e.g. in the diagrams of 
the {4, 2} type. However, before the final integration over 
the internal coordinates, we seek contributions, which are 
proportional to the first power of external frequency f2„, 
this linear dependence coming from two external vertices 
belonging to the same loop. If these vertices belong to 
two different loops, then such diagram is proportional 
at least to il^^. After analytic continuation iil„ — > 
and (non-logarithmic) integration over interaction region 
(—L,L) we have an extra factor, HL, which is small in 
the limit — > 0. In other words, we neglect the screen- 
ing of external field, which is given by RPA series where 
fermion polarization loops are (1-reducibly) connected by 
interaction lines. It is allowed when measurement times, 
^ , are large compared to times of travel through the 
interacting region, L/vp. Thus we find the value unity 
for the conductance in the absence of the barrier i^iii 



B. From perturbative corrections to the RG 
equation 

In this subsection we discuss the general form of cor- 
rections to the conductance and outline basic ideas of 
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the renormalization group approach, which we will use 
below. 

Starting with the conductance < G < 1 in the non- 
interacting limit, the corrections 5G induced by the in- 
teraction take the form of a power series in g. More 
precisely, we discuss the equivalent quantity F = 2G — 1 
, equal to cos26' in the limit g ^ Q . We denote the 
quantity Y , renormalized by interaction effects, by 
and recall that it may be represented as a power series 
in g . The coefficients of this series contain terms vary- 
ing as powers of A = ln(L/Lo) (at zero temperature) 
or A = ln(To/r) at finite temperature T in the limit of 
large L . Here Lq and Tq are ultraviolet cutoff param- 
eters. Accordingly, Yr may be represented as a double 
series. 



34)) with respect to A and put A = 



Yr = Yo + biigA 

+b33g^A^ + 1)3.29^ 



(32) 



with hij functions of Y 



The terms of the form (.gA)" 
in l|32p are conventionally called "leading" contributions, 
and the terms of the form 5" A™ are called "subleading" 
for n > m. The first term Yq is the sum of all scale 
independent contributions 



Fo = i;(A - 0) - r + 52o.g' + 630.9' + 



(33) 



Naively speaking, the correction should be small SY = 
Yr — Y ^ Y, in order to be considered as a correc- 
tion. More precisely, the series ((32)) should be converg- 



ing. When all \bij\ ~ 1 and A 1, this convergence is 
achieved for gA <C 1 which in turn imposes stricter con- 
ditions on the smallness of g. However, assuming that 
Yr{g,A) is an analytic function of A, we expect that an 
analytical continuation from the region of small A to large 
A should be possible^'^ . Then we can relax the require- 
ment of the smallness of g. Rearranging the above series 
we write 



y,(A)=yo+/3iA+yA' + §A' 



(34) 



where the coefffcients (3i are functions of Y . 

If the theory is renormalizable, the following scaling 
relation holds: 

Yr{g,Y,T) = F{g,T/e) 

where Q = 8(5, y) is the correlation temperature of the 
problem (an analogous relation holds for the scaling with 
L at zero temperature). The scaling property implies the 
existence of a renormalization group equation 



dYrig,Y,A) 
dA 



^P{g,Yr{A)) 



(35) 



where f3{g,Yr) is the so-called /3-function. To obtain f3 
from perturbation theory we take the first derivative of 



dYrig,A) 



dA 



Pi{g,Y) ^ P{g,Yo) 



(36) 



A=0 



In the last equation we have replaced Y by Y{Yq) , ob- 
tained by inverting the series ()33p for Kg in powers of 
y, which transforms the function /3i(g, Y) into (3{g, Yq) . 
Now we see that (|36)) is nothing else but the RG equation 
(|55)) at A = where Yr — Yq. Consequently, the function 
/?((7,yb) is the true /3-function. 

To check whether the scaling and therefore the RG 
equation holds, we may integrate (|35)) to recover the per- 
turbation series, by employing the inverse function 



A = 



dY 



(37) 



Alternatively, noting that ^ = f3{Y)-^ we may obtain 
the Taylor expansion in powers of A directly in the form 



YriA) = Y + Ap{Y) + -A^f3iY)p'iY) 

+ iA3/3(r)[/3(y)/3'(r)]' + ... 
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(38) 



with a prime denoting the differentiation with respect to 
Y. 

In what follows we, firstly, verify the structure (|38)) of 
the theory up to third order in g, using computer algebra. 
Thus we show the applicability of the RG approach. Sec- 
ondly, on the basis of this analysis we propose a method 
of resummation of the most important contributions to 
the /^-function in all orders of g. 

But before presenting these results we return to the 
significance of the scale independent terms in the pertur- 
bation theory (the difference between Yq and Y), as this 
problem has been sometimes overlooked in the literature. 
In the first order of g we can absorb such finite terms into 
the definition of A. However, after this fixing of the defi- 
nition of A, higher order scale independent terms cannot 
be removed by redefinition. For example, we should ex- 
pect a term 621^^ (A -I- c) with c ~ 1 in the second line of 
(|32)) . It is clear that in the third order we can find a corre- 
sponding multiplicative contribution biib2ig^A{A + c), 
which contains a linear-log term and thus affects the def- 
inition of the /3— function in the form of (|36)) . Evidently, 
such terms arise only at the level of g^, or "in the third 
loop" , and do not affect the leading- log contributions. As 
this problem is of principal importance, we present in the 
following a second line of derivation of the /9-function. 

Within the Callan-Symanzik (CS) formulation of the 
RG metho d^^'^^ one may start with the perturbation se- 
ries of Yr in terms of the bare conductance Y , ((32)) . The 
latter relation is inverted and Y is expressed in terms of 
the "running" parameter Yr 



Y = Yr + bngA 

+b229^A^ + hig^A 
+b33g^A^ + &325'A2 



&31.9 A 



(39) 



10 



with hij functions of Y^.. The expansion ((39)) is the essence 
of the CS formulation, appropriate for our problem. Here 
Y is understood to be a function of the two variables 
Yr{K) and A . Since the quantity Y is independent of A 
we have the relation 



= 



dY_ 
dA 



dY 



dY dYr 
dYr dA 



or 



f3 



dYr _ dY/dA 
'dA ~ " dY/dYr 



(40) 



(41) 



In practice, the calculation of the /3-function from ([1T|) 
requires the calculation of the coefficients bij from pertur- 
bation theory, which is analogous to the transformation 
of Pi{Y) to /3{Yo) discussed above. If everything is done 
correctly, then the expression (|4T|) does not depend on A 
explicitly. 

The CS procedure described above is different from 
the Gell-Mann-Low (GL) formulation of RG. In the lat- 
ter one considers the value of conductance used in the 
calculation to be defined at the scale A ; which means Y^ 
in our notation. To eliminate the large terms involving 
A one introduces counterterms into the Hamiltonian (cf. 
below), in order to guarantee a finite value of the bare 
conductance, with a resulting equation similar to ([M)) . 
Usually the insertion of the counterterms is done in the 
"minimal subtraction" scheme, meaning that scale inde- 
pendent contributions are omitted. As we show in Ap- 
pendix [Gl this GL formulation leads to ambiguities in 
the definition of the /3-function at the level of the third 
loop, in contrast to the CS formulation or the scaling 
formulation presented in Eq. p6p . 



C. Summary for the zero temperature case 

In the perturbation theory calculation the logarithmic 
factors appear at the last step, when integrating over the 
points of interaction. Prior to this step the calculations 
are straightforward. We proceed first with the zero tem- 
perature case. In this case the sums over discrete Mat- 
subara frequencies transform into integrals over imagi- 
nary energies, which are rather easily done by computer 
symbolic computation. 

Each diagram contains a prefactor $76^^^'^^+^' with the 
two external vertices at x,y > 0. Integrating over y and 
putting 17 — > makes this prefactor unity. Equivalently, 
we may keep only the first, linear-in-r2, term in each dia- 
gram and let r2 = in its remaining part. In the first or- 
der of g the correction is given by the diagrams in Fig. [2l 
The vertex correction part in Fig.O; cancels and the two 
first diagrams give an expression which contains 1/xi. 
This should be integrated over the point of interaction 
Xi in the limits (a, L) which produces a factor 






FIG. 2: Three Feynman diagrams, depicting the first order in 
g2 contribution to the conductance. The other three diagrams 
are obtained from these by reversing the direction of fermionic 
propagation. 



Similarly, in the second order of g we have two contribu- 
tions, {xiX2)~^ and {xi + ^2)"^, which lead to Ag and 
Aq, respectively. We provide further details in Appendix 
IeI and list here only the summary of our calculation. 

Grouping corrections in powers oi g, as in (|32|) we have 
for the renormalized quantity Yr'. 



Yr 



= y-C{4}5A(l-F') 

+ g^Y{l - Y^)[-C[e}{A^ - 04) -f C{4,2}(A - a,)/2] 
+ 5^(1 -Y^)[~ C{8}(A3 - 3Aa4)(>'' - 1/3) 
+ C{8}(l-r2)Aa2+C{6,2}i^'A(A-ai) 

- C{4,4}(l-i"')(AV4-a3A) 

- q4,2,2}(l + Y')A/4 



(43) 



where a = (21n2, 7rVl2, (ln2 - l)/2,0). AU c 



{■} 



1 



and these coefficients are shown only for reference to the 
underlying group of diagrams, according to the above 
conventions. 

Using the above CS scheme we come after some algebra 
to the RG equation dY^/dA = (3{Yr) with 



[l-Y^) 
C35'(l - 



>Y 



>Y^ 



1 
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Y 



2\2 



2 

0{g') 

-2 l-ln2 



(44) 



C3 = 02 



03 



04 = 

12 



Aq = ln_L/a 



(42) 



Here the value of C3 is that of the T — case (cf. below) . 
This value differs from our previously reported resuHi^l, 
which accounted only for the set {8} of third order dia- 
grams. Here we found an aditional contribution from the 
{4,4} set. 

It is important to keep the scale independent contribu- 
tions in the terms {6} and {4, 2} in in order to have 
the agreement with the result by Kane and Fisher in the 
limits of weak scattering and weak tunneling Y"^ 1. 

One can check that the terms with higher powers of A 
in (|43|l are reproduced by Eqs. ([38]) and (j44|) . In fact, this 
applicability of the RG equation is manifested already 
in the absence of A in the right-hand side of (|^^ . This 
means that, when discussing diagrammatic contributions 
in nth order in g we should concentrate only on the least- 
leading logarithm, i.e. the terms ~ g^A. These linear-in- 
A contributions arise in two qualitatively different cases. 
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+ symm. 



FIG. 3: Skeleton Feynman graphs, leading to lowest order 
logarithmic contribution, g^A, beyond ladder series. 



In the first case, the entire Feynman diagram is pro- 
portional to a single logarithm. This happens with the 
three first terms in (j44p stemming from diagrams with 
the maximum number of fermionic loops. In the above 
notation, these are the sets {4}, {4,2}, {4,2,2} for cor- 
rections in order g, g^, g^, respectively. We discuss these 
corrections in the next section. The finite terms, 0(1), 
do not contribute to the /3-function as can be checked by 
the absence of coefficient ai in p4)) . 

In addition to these sets, there are contributions in the 
third order, whose leading divergence is also linear loga- 
rithmic. They are depicted by a first skeleton graph in 
Fig. [21 with the maximally crossed wavy lines of interac- 
tion. These graphs provide for one half of the value of 02 
in Eqs. (031), (glD. 

In the second case, the linear-log contributions arise 
as accompanying weaker divergence in the stronger di- 
vergent graphs. In simple cases, it may be a combina- 
tion -I- AO(l) ; in more complicated situations, it 
may be a stronger singularity at internal points, e.g. 
~ J dzi/{zi — Z2)~^ which is eventually removed upon 
symmetrization, whereas the subleading log-divergence 
survives and contributes to ([321), cf.^®. This second 
type of linear- log contributions happens in the remaining 
graphs in Fig. [3l and provides for the rest of the value of 
the last term C3 in 



D. Finite temperatures 

Let us discuss the effects of temperature, first in the 
lowest order m g. At T = the large logarithm of the 
theory (|42p came from the integration J dxi/xi over the 
range of interaction [a, L\. When the integration over in- 
ternal frequencies is replaced by a summation over Mat- 
subara frequencies, we find that the integrand contains 
the fermionic Green's function Q{2x,t = 0). of the form 



T 1 
2vf smh{2T:Tx/vF) A'kx 

with an appropriate change in the definition 



An 



.'^ , X , tanhL/£ 

47r / dx g{2x) = In - — — ^ = A 



tanh a/£ 



(45) 



(46) 



with the temperature correlation length, £ = vp/nT. Ex- 
plicitly we have 



A = lncoth(a/£) 



(47) 



where the last relation takes place at larger temperatures, 
when L ^ £. Here we introduced the ultra-violet energy 
cutoff To — vf/t^o-; one may think of a ~ kp^ so that 
To ~ Ep- The results presented below refer to this case 
of relatively high temperatures, Tq^ T :§> vp/L- 

We calculated the terms of perturbation theory by 
means of computer algebra, similarly to the T = case 
above. The summary of the result is given by the same 
expression ([151) . but the subleading coefficients are now 
given by a = (3 In 2, TT 2/12, In 2 - 1/2, 7r2/24) . In its turn, 
it leads to a different numerical value of the coefficient 
([ill) , we obtain for T > uj^/L : 



C3 



tt' 1 
24+1-2-2 



(48) 



Let us comment on this result. 

We observe that the diagrams with the maximum num- 
ber of fermionic loops lead to expressions linear in A with 
the same prefactors as in the T = case. They generate 
the same part of the /3 function, given by the first line in 
(|44p . and correspond to one-loop RG approximation. In 
the next section we show that these terms form a series 
which can be resummed in all orders of g. 

The maximally crossed diagrams of the third order 
which had only linear-log divergence at T = lead to 
the same 02/2 coefficient in ([151) a.t T ^ vp/ L. 

Let us now consider the modifications in the second 
order. The complete expression for the diagrammatic set 
{6} is given by the expression (|F2p which shows that 04 
in Eq. (|43p is a smooth finite function of temperature. 

At finite T each particular diagram in order g^ can lead 
to a rather complex expression, and the resulting sum of 
all diagrams becomes too complicated for a brute-force 
approach. The analysis can be best done in the follow- 
ing way. We add all the terms of a given structural set 
of diagrams, thus removing the internal singularities at 
zi = Z2, etc. As discussed above the leading logarithmic 
behavior of the set (higher power of A) is determined by 
the RG equation. We may hence subtract a simpler ex- 
pression, with the same leading behavior, from the com- 
plicated general expression. 

For instance, for the {8} set we pick all di- 
agrams proportional to cos 80, and notice that 
the zero-T expression has a leading contribution 
(21222:3)^^ with a certain coefficient. Then we ex- 
pect that the finite-T expression should have a leading 
term 8C"3(sinh(2zi/£) sinh(222/C) sinh(2z3/C))"^ with 
the same coefficient. We subtract this term from our 
expression and analyze the degree of singularity of the 
remaining expression. It turns out that in the third or- 
der after the subtraction of such leading A"^ terms we are 
left with linear- log expressions. The coefficients before 
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these expressions are also smooth functions of T, as we 
explain in the Appendix [Fl 

Concluding this section, we verified the validity of 
the RG-equation by checking that the higher-than-linear 
powers of A are generated by the linear-log terms. This 
holds both at r = and at T ^ vp/L; the resulting 
/?-function is somewhat different in these cases and we 
discuss this difference below. There is a part of the (3- 
function, which is identical in both cases, and stems par- 
ticularly from a sequence of diagrams with the maximum 
number of polarization loops. We analyze this sequence 
in the next Section and show that it is possible to resum 
it in all orders of g. 

It is important to note that the last term in ((44|) is van- 
ishing in the limit when the barrier becomes fully trans- 
parent or fully reflecting, i.e. at y ^ ±1. In order to see 
that, let us consider the case (y+l)/2 = G ^ 0, then the 
first three terms in (gi]) read (3 ~ 2G [-g - g^/2 - g^/2] 
and lead to power-law scaling of G in accordance to re- 
sutls by Kane and Fisher (see below). At the same time 
the last term in (|^^ is oc g^G^ and does not lead to a 
noticeable effect in the limit considered. 

It is also worth noting here, that the calculation both 
at T = and at T ^ vp/L confirms the relation between 
the subleading coefficients in sets {8} and {6,2} on one 
hand and {6} and {4, 2} on the other hand. This is seen 
in Eq. (j43|l for the a4 and ai coefficients. Without this 
relation the results by Kane and Fisher for weak and 
strong barrier would not be recovered. 



VI. COUNTERTERMS IN THE HAMILTONIAN 

In the previous section we saw that the corrections to 
the conductance are obtained as self-energy insertion into 
the fermionic loop, whereas the vertex corrections vanish 
in the dc limit Q <^ vp/L. This means, in particular, 
that the effect of renormalization in this limit can be 
found at the level of the one-particle Green's function, or 
by addition of counterterms to the effective low-energy 
Hamiltonian. 

We calculate the corrections to the matrix Green's 
function for the chiral fermions in the pseudospin sector 
in perturbation theory. The rules of the diagrammatic 
technique remain the same. This time we do not con- 
sider a closed fermionic loop, which described the con- 
ductance above. Instead we study the matrix Green's 
function G{x, —y, ri„) connecting distant points on differ- 
ent sides of the barrier in the chiral representation. The 
initial Green's function is diagonal in pseudospin space, 
and the corrections to it make it non-diagonal. We show 
that this non-diagonal form corresponds to additional ro- 
tation of pseudospin near the origin, in accordance with 
consideration of Sec. IIVI 

Employing the CS scheme we start with the bare bar- 
rier, characterized by the parameter 9 and find the renor- 
malized value 6', (0, A) as described below. Then we de- 
termine the inverse relation 0(0^, A) and require the in- 



dependence of initial 6 at the scale of consideration A. 
This provides us with the RG equation for 6r 

The prefactor in G{x, —y, fin) — G exp — il„(a; + y) is 
found as 



G = 1 + C|4}^io-i sin20 

-C{6}^(sin2 26* A^ - iai sin46'(A2 - bi)) 



(49) 



40 



-C{4,2} sin46'(A - 63) 

3 

-C[s} 1^ (3(A^ - 61 A) sin 20 sin 
A3 

-icTi— sin26'(9cos46l- 1) 
+il2cri A sin 26'(&2 sin^ 26* + 61 ) 
+C{6,2} A( A - 63) sin 46'(sin 26 - 2iai cos 26) 

+Cf4,i}^icri sin^ 26'A(A - 64) 
8 

+C{4,2,2} ^^o"! A(5 sin 26 + sin 66) 

where again all C{.} = 1. The coefficients = 1,...4 
depend on the case in consideration. For the zero tem- 
perature we find h = (0, 7r^/6, 2 ln2, 2 — 2 ln2) and for the 
finite temperature 6= (7rVl2, 0, 3 ln2, 2 - 41n 2). 

In order to establish the connection to the Hamiltonian 
we observe that the above form ([49]) can be represented 
as = exp ia\ {9r — 6) with 



0r = e + ^sm2e + ^smA0{A^ -A-bi + bs) 
2 8 

-h— sin26'f-A3cos46i + A^{1 + 3 cos 46*) 
16 V3 

-hA((3- 461-262 + 263 -64) 

+ (1 + 262 + 263 + 64) cos 46*)) + 0{g^A°) (50) 

This means that the effect of the interaction is indeed re- 
duced to a change in the magnetic field 20, in accordance 
with Sec. [mP] 

We have for the value at A = 

^'o = 6' + ^(&3-fci)sin46l + ..., 
8 

Expressing by 0o in (jSOp and inverting the series we 
obtain 

6*0 = Or-^ sin 261^ -I- ^ sin 4:0r{A'^ + A) 
2 8 



16 



■ sin 20r (- A^ cos 46*^ + A^ (1 3 cos 40^) 
+A{{3 - 2bi - 262 - 64) 

-|-(l + 26i-^ 262 -^64)cos46'^)W... (51) 
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The difference Oq — Or: expressed in terms of Or, is inter- 
preted as the counterterms to the Hamiltonian, needed 
to compensate the divergent diagrammatic contributions. 
In our approach we do not re-calculate diagrams with in- 
clusion of the appropriately chosen counterterms, which 
would be problematic, given the non-Abelian character 
of the theory and the absence of Wick's theorem. These 
counterterms merely reflect the structure of the correc- 
tions to the magnetic field. Sec. IIIIDI above. 

Demanding the independence of on A, we obtain 
(here 6'(A) = Br) 

89 0^ 
2— = gsin26'- — sin46'-|- — sin26'(l -|-cos2 26') 
OA 4 4 ^ ' 

sin3 20 0(g4) (52) 

where 



loops. These terms are cx (1 — Y'^) and therefore are 
important at any value of the conductance. It is possible 

-X\ -Xi -Xi -X, -Xi -z -X2 -Xi -z -X2 



X2 




Xi Xi 



FIG. 4: Ladder series, L{xi,X2;ojn), describing the combined 
effect of interaction and barrier, and leading to a linear-in- 
logarithm ln(ro/r) contribution to conductance. 



to resum all the contributions of this type, as we discuss 
below. 



= -(261 + 252 + 64) 
7r2 ln2 1 



T = 



24 



ln2- 



1 



0.60, T:$>vf/L 



which agrees with the above Eqs. (j44|l . (j48|l . This coin- 
cidence confirms that the vertex corrections to the con- 
ductance are unimportant in the considered dc limit. 

Any method of rendering the logarithmically divergent 
integrals finite is called a scheme of regularization. As we 
explained above the contributions linear in A get addi- 
tional contributions in third order of g from the scale 
independent terms appearing in the second order of g. 
We removed this ambiguity by using the scaling method 
or the appropriately devised CS scheme. A part of the 
answer is independent of the regularization scheme and 
technically is obtained by summing the one-loop contri- 
butions to /5-function. One loop means here one inde- 
pendent integration over the internal energy, leading to 
a linear logarithm. Indeed the summation of the ladder 
series in Sec. I VIII amounts to a finite renormalization of 
the effective interaction constant, g g- In addition to 
that, the quantities calculated above in the order g^ in- 
cluded linear-log contributions stemming from diagrams 
with three independent energy integrations (Matsubara 
frequency summations). Particularly, the set of diagrams 
{4, 4} contains diagrams with vertex corrections to the 
polarization loops in the ladder series Fig. 21 which con- 
tribute to the rcgularization-dependent coefficient 64 in 
Eq. (gni) and C3 in ^ 



VII. 



INFINITE RESUMMATION OF TERMS IN 
THE /3-FUNCTION 



A. Ladder series and Wiener-Hopf equation 

The above discussed linear-log contributions to the 
conductance with the maximum number of fermionic 
loops, can be viewed as a "dressing" of the wavy interac- 
tion line g in Fig. [5^, Fig. by fermionic polarization 
loops. This is shown schematically in Fig. 01 with the 
result of the summation denoted by g and depicted by a 
double wavy line. 

We show below that the result of this summation is 
scale independent (without logarithmic terms), so that 
using the double wavy line in first-order contributions 
Fig. [5^, Fig. [2b will produce a single logarithmic factor, 
while containing all powers of g. 

In the usual situation, the RPA summation of Fig. 3] 
is trivial after Fourier transform. In our case of chiral 
fermion representation, we have left-right asymmetry of 
the vertices whose value in pseudo-spin space depends 
on the coordinate. Hence we always integrate over the 
semi- axis only, which leads, as shown below, to an in- 
tegral equation of the Wiener-Hopf type, instead of a 
simple Fourier convolution. Still, the ladder equation for 
the dressed interaction can be written and solved rather 
simply. 

Let us denote the interaction between the points —x 
and y by L(x, y; ujn) where tUn is the Matsubara frequency 
along the line. Let us define the vector X{x, y. uJn) in the 
following way 



X{x,y;u„ 



_L{x,y;uJn) 
L{x,-y;ujn. 



(53) 



We observe that the first three terms in (|33]), ([5^ are 
given by diagrams with maximum number of fermionic 



Then the ladder summation. Fig. [U is expressed for 
> by the integral equation : 
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X{x,y;un) 



-■iirgSix - y) 



dz 



.X{z,y-ujn) (54) 



where we used the fact that the polarization loop in the (w„,a;) representation is {J2,{x,T)J'i{y,Q))^^ = 
(47r)~^a;„0(aj„cc)e~'^"^. Here and below the integration over coordinates is over the interval (a,L), it is convergent 
and may be extended to (0, oo). 

Later we will need the integrated quantity L{x;uJn) ~ J dye^'^"^ L{x,y;LOn), which obeys the integral equation 



47r + w(y + |) / dze-'^^L{z]Uj) 



+ ^g^uj J dz e-"l^-^lL(z;w). 



(55) 



B. Solution and its properties 

The above equation for L{x; uj) is an integral equation of Wiener-Hopf type, and its solution can be found as follows. 
We differentiate both sides of Eq. ([55]) twice with respect to x and find that d^ L{x] to) / dx^ = a;^(l — g^)L{x;uj). 
Demanding that L{x ~* oo;lu) = we seek a solution in the form L(x;u>) — C exp(—ujxy^l — g"^). Substituting the 
latter form into (I55p . we determine the constant C and find 



L(x; uj) — —4:Trg 



1 + VI 



1 + Vi^ 



(56) 



Now we are in a position to determine the linear-log correction to the conductance stemming from the diagrams Fig. 
[2] where the wavy line of interaction is replaced by the double wavy line. 
The equation reads : 



dxidx2dy L{xi, X2]uj)Q{xi - x; e)Q{-xi - X2;e- iLj)Q{x2 - y; e)Q{y + x;e + fl) , (57) 



where the superscript L in G^^^ denotes ladder summa- 
tion. Taking here the limit 17 ^ we find 



above, we find 



^(i) ^ -2ff(i-y^) 



l + ^l-g^+gY 



Ao = -5(1 - Y^)Ao . (58) 



Here we defined the renormalized interaction constant g. 
We observe that 



L{x,y;ujn) 



-5{x - y) 



2d 



-ujnd\x — y\ 



(60) 



-u,.d{x+y) Y{l+d)+g 
2d Y{\ -d)+g 



1 - K, 



Y : 

Y 



(59) 



with the Luttinger parameter K — — g)/{l + g). 

In fact, the expression (l57|) already assumes the limit 
17 ^ 0, because the prefactor 1 ~ Y^ is obtained by 
adding two contributions, one shown in Fig. [2j3 (cx $7(1 — 
2Y'^)) and another in Fig. [2^ with inverted direction of 
fermionic lines (oc 12). 

The importance of the limit 17 — > in the derivation 
of ([58]) is further illustrated by the solution Z(x, y; w„) of 
the initial Eq. ([5^ rather than Eq. . After a long and 
straightforward calculation, similar to the one described 



with d = yjl — g"^. The first term in (jHO]) is the bare 
interaction, the second one is the result of RPA summa- 
tion, which is independent of the barrier transparency Y 
and finite even in the bulk, at x ~ y ^ oo. The last term 
explicitly contains the transparency Y and decays away 
from the barrier. Integrating over y we get the above 
Eq. ([5S)) . and Eq. ([Sn]) shows that the "dressing" of the 
interaction constant, g — > g, in (j58p is not reduced to a 
simple multiplicative factor, but depends on the coordi- 
nates. Finally we note that G^^' = at A = 0, i.e. there 
are no scale independent contributions in the ladder ap- 
proximation, and Yq —Y \n that case. 
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VIII. RG EQUATION AND ITS SOLUTION 



Integrating Eq. ([65)) we obtain the following result : 



From the solution of the Wiener- Hopf equation for the 
ladder series (|55| we then find the /3-function in ladder 
approximation 



Pl{Y) = - 



2g(l-y2) 



1 + yi" 



(61) 



Inverting this equation and using the definition of K we 
obtain the symmetric form 



dA 

dY 



1 



I- K\-Y 1 - A'-i 1 



Y 
(62) 

invariant under the simultaneous exchange oi K K ^ 
and y <-> — y . This differential equation may be solved 
analytically-^^. We recall that it is asymptotically exact 
in the limit |y| ^ 1 . 

Let us now turn to the first correction beyond this 
series in Eq. ([33]). Applying the correction as calculated 
we get the improved RG equation 



dY _ 2^(1 

rfA " 1 + ^1 - g2 + gY 



+ C3.g3(l - y2)2 (63) 



The correction term is qualitatively different in that it 
vanishes quadratically in the limit |y| — > 1 and hence 
is irrelevant in this limit (considered in the early work 
by Kane and Fisher). The duality property still holds, 
i.e. the invariance of (|63p under the simultaneous change 
g ^ -g and r ~Y (or e 7r/2 - 6* in Based 
on these observations, we expect that the fourth-order 
terms beyond the ladder series in the RG equation will 
be proportional to g'^Y{l — Y'^Y . One should keep in 
mind that the correction terms here have been included 
only to order . We now discuss a further improvement. 

We may perform yet one more partial resummation by 
dressing interaction lines as g ^ where 



l + K 
1 - K 



Y 



(64) 



The C3 term in (j63p results from non-ladder diagrams. 
Therefore without double-counting we can dress the in- 
teraction constant there and write 035^(1 — Y'^Y . Then 
([S^ receives a correction term of the form 



dA 

dY 



1 



1 



C3 



1 - a: 1 - y 

2(1 -i^) 

1 + A' + y(i - A') 



1 - a:-i i + y 



(65) 



The term proportional to C3 here is of order of g at small 
g and we omitted higher order terms. 



$(G)/<i>(Go) 



1-G 



(66) 



we assume the initial condition G = Go = (l + yo)/2 = 
cos2(e') -I- C3gocos(26')sin2(26') at T = Tq, and g^^ ^ 

i i^jz^ + cos(20)^ ; here C3 is given by ([i5|l for the case 

T 3> vf/ L. Evidently, the renormalization stops at T < 
vf/L. 



IX. SUMMARY AND DISCUSSION 

In this work we have presented a fermionic formula- 
tion of the transport theory of interacting electrons in 
a onedimensional system with potential barrier in the 
linear response regime. We considered the Luttinger liq- 
uid model with interaction 32 (forward scattering involv- 
ing right and left movers), and employed a current al- 
gebra represention, leading to a substantial simplifica- 
tion of perturbation theoretical calculations. Using the 
renormalization group idea in various forms we showed 
that the conductance as a function of temperature (or, at 
zero temperature, as a function of system length) obeys 
scaling, and therefore may be obtained by integrating a 
renormalization group equation. The /3-function of the 
latter equation has been found within a ladder approx- 
imation, which becomes exact in the limit of small or 
large conductance, at any interaction strength. Correc- 
tions to the ladder approximation up to and including 
third order in the interaction, and relevant at interme- 
diate values of the conductance, have been calculated. 
These correction terms require the solution of a funda- 
mental problem of renormalization group theory: What 
is the significance of and how does one treat scale inde- 
pendent contributions appearing in perturbation theory? 
We show that these terms have to be taken into account, 
as they guarantee the scaling of the conductance and can 
lead to an important redefinition of the /^-function. In 
that respect our work is of general interest in the con- 
text of the renormalization group method, beyond the 
particular transport problem considered here. 

Eq. ([5S)) is one of the central results of our paper. It 
is in agreement with the previous findings^i^ii^ in the 
limiting cases G ^ 1, G ^ (i.e. |y| ^ 1 ), AT 
1, except for the fact that in these previous works the 
conductance tends to K, rather than 1 , in the limit of 
vanishing barrier strength at finite temperature. In this 
limit ( |y I ^ 1 ) our result is exact for any value of AT . In 
the intermediate regime G ~ ^ there appear corrections 
to our result. We have calculated these corrections in 
the lowest order, which is third order in the interaction 
strength. The parameter C3 determines the strength of 
these corrections. We have shown that C3 is nonuniversal 
in the sense that it depends on the physical nature of 



the ultraviolet cutoff, e.g. the temperature cutoff or the 
length cutoff. We are convinced that the scaling functions 
of the temperature dependent conductance, G{T), in the 
limit of infinite sample length L and the length dependent 
conductance, G{L), at T = are slightly different. 

Let us discuss the role of the factor C3 in (|66|1 which 
is new as compared to previous studies^>^ii^ . It is clear 
that C3 in ([5^ defines the renormalization of the con- 
ductance in the intermediate region, G ^ 1/2, and thus 
determines the coefficient relating high-temperature and 
low-temperature asymptotes. We may ask the following 
question: what is the coefficient C in the asymptotic ex- 
pression G ~ Cr2(^"'-i) (r = T/To ) if we fix the overall 
temperature scale at hight T as G ~ 1 — r^i^-^) ? From 
Eq. (inH) we have 



G = ii:4c3(i-i/-ff) 



(67) 



At this point it is appropriate to compare our findings 
to the exact results obtained within the conformal field 
theory (CFT) approach to the bosonized version of the 
problem. Sec. Ill Al In these earlier works the /3-function 
of the RG equation was not discussed. However results 
equivalent to those by Fendley et alJ^^ are presented in 
a recent study by Lukyanov and Werner—, where partic- 
ularly the /3-function was obtained in a form equivalent 
to our (I66|) . According to Lukyanov (private commu- 
nication), the coefficient C3 = 1/4, which deviates from 
our value C3 ~ 0.6 and we use both values for compar- 
ison below. Unfortunately, one rarely, if ever, discusses 
the scheme of regularization used in obtaining the CFT 
exact solutions. We have to conclude that the regular- 
ization used bji2^ is different from ours. The choice of 
cutoff scheme in our work is not based on mathemati- 
cal convenience, but rather it is following naturally from 
perturbation theory. Therefore, as a theoretical result 
to be applied to explain experimental data, our (|66[) has 
advantages over the CFT results. 

We saw above that passing from the backscattering 
amplitude in the Hamiltonian to the conductance (i.e. 
5-matrix) may require ultraviolet regularization, i.e. cer- 
tain prescriptions of dealing with short-distance singular 
quantities. Our observation is supported by the study 
of the free {K = 1) BSG theory^i^ where the authors 
arrive at a current algebra formulation similar to our 
treatment above. They notice that the description of 
the strong coupling limit \u2\ — > C!0 corresponds to tak- 
ing the quantity \9\ — > tt, (see Eq. (2.23) there), which is 
related to our |B| = 29 and hence to the conductance. 
They also notice that the connection between \9\ and 
\u2\ = 9f{9^) ^ 9 + 9^0{1) begins to depend on the 
type of regularization chosen in the third order of 9. On 
the other hand, given that the CFT exact solution in 
the interacting case if ^ 1 is expressed in terms of the 
scaling variable |u2|(7o/7')^~^, we see that the implied 
equation \u2\{To/T)^~^ ~ 9f{9'^) corresponds in form 
exactly to our (|66p. It would be thus desirable to un- 
derstand the particular regularization used in the CFT 
solution in more detail. 
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FIG. 5: (Color online) (a) Conductance as a function of tem- 
perature for K — 1/2. The results of ladder summation, (|66p 
with C3 = 0, are shown by dashed-dotted line ; the result of 
Eq. (|66p . with C3 given by ((52} is shown by solid line. The 
result available from CFT (multiplied by factor 2) is shown 
by pluses.— i^ii^, together with our fit by (|66p with C3 = 1/4 
shown as dashed line, (b) The same plot in log-linear scale. 



Let us now illustrate the details of the behavior of G(T) 
in the example of the well-known case K = 1/2. Both 
high-temperature and low-temperature asymptotes of G 
are interesting from the theoretical viewpoint, but one 
can expect the experimental variation of the conductance 
to be most visible at lower T, when G < 1/2. It is there- 
fore useful to consider the variation of the conductance 
for fixed K and different values of C3. We take K = 1/2 
and and fix the low T behavior in the form G :^ , with 
T = T/n and ~ To|r|2. We then plot G(r), obtained 
from ([M]) for several values of C3 in Fig. [51 together with 
the exact solution, known froniii^iiS 



Gl/2 — 1 - 



— ^ip' - H ^ 

2V3r V2 2\/3r 



(68) 



where tp'{x) the derivative of digamma function. For K = 
1/2, and the implied CFT value C3 = 1/4 we can invert 
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(l66l) and write 



G 



2r2 



1 + 2r2 + VI + 8t2 



(69) 



The Fig. O shows that the CFT solution Ues between 
the curve corresponding to pure ladder summation, C3 
0, and the one with the above C3 ~ 0.60. The proposed 
value C3 = 1/4 fits the CFT solution rather precisely. 
This tells us that the approximation leading to (j66p works 
rather well at the interaction strength K — l/2{g — 3/5) 
, provided the coefficient C3 appropriate for the regular- 
ization scheme used in the CFT solution is taken. In fact, 
we show in the Appendix |H1 that the overall agreement 
criterion using the coefficient C, Eq. ([57]) . is about 1 % in 
this case. For a further case reported {K = 1/3 , g = 4/5 
the agreement is worse, about 10 %, which might 
be expected in this more strongly interacting situation. 

Finally, we emphasize again that in any comparison 
with experiment the correct regularization scheme is de- 
termined by the details of the corresponding experimen- 
tal setup. Indications of the possible inadequacy of the 
regularization scheme used in the CFT solutions exist. 
For instance, we tend to interpret the systematic devi- 
ation of the experimental data on the conductance be- 
tween the quantum Hall edge states from the CFT pre- 
diction at low T in^-. not as a shortcoming of the ex- 
periment, but rather as a consequence of choosing an 
inappropriate regularization. The correct regularization 
might lead to a different value of C3 in (|66p . providing 
better agreement with experiment. 



Acknowledgments 

We are grateful to I.V. Gornyi, D.G. Polyakov, P.M. 
Ostrovsky, K.A. Matveev, D.A. Bagrets, CM. Yev- 
tushenko, M.N. Kiselev, A.M. Finkel'stein, A. A. Ners- 
esyan, L.I. Glazman and N. Andrei for various useful 
discussions. 



APPENDIX A: FROM THE HAMILTONIAN TO 
THE S-MATRIX 

The continuum Hamiltonian ([3]) corresponds to a lat- 
tice model 

00 ^2 
Hiat = X! + h.c.) + 2 + (-l)^"2)c]cj, 

at half-filling (2fci? — tt) and with potential present at 
two adjacent sites. 

In Eq. ([3]) we let ui^2{x) — wi,2<5o(i,2)(a;)- Here ui,2 
are dimensionless quantities and 01^2 is the range of the 
corresponding regularization of the 5-function; at the end 
the limit oi 2 ^ will be taken. 



Combining the two chiral fermions into a spinor quan- 
tity = (V'jf, V'l) we can write the equation of motion 

as 



\ U2{x) -idx+ui{x)J ' ^ ' 
VF{i<7-idx + ui{x) + U2{x)ai + (x)cr2)^', 



with real and imaginary parts U2 = u'2 + 11*2'. The formal 
solution of (|A1|) is 

*(x,t) = e*"*r(a;,a)*(a,0), (A2) 
T{x,a) = %expi / dyaslvp^uj 

J a 

+ui{y) + u'2{y)ai + 4'(y)a2], (A3) 

where stands for the x-ordering operator. In the long 
wavelength limit, w ^ 0, we obtain 

T{x,a) = exp[iuJVp^a3{x — a)], x<0, (A-4) 
= exp[iujVp^a3x]To exp[—iujVp^a3a], x > 0, 

with To depending on ui^2{x). We assume first that the 
regularized J-functions in the definitons of ^1^2(2^) are 
identical, which means ai =02- In this case 



To = exp(iui(T3— M2CT2-f U20'i) — coshb+bcr 



sinhfo 



(A5) 



with b = (7/2, —u'2, iui) and — |M2p — u\. This trans- 
fer matrix simplifies in two important cases : i) point-like 
impurity, when ui — \u2\ and ii) purely backward scat- 
tering impurity, ui — 0, U2 ^ 0. We have 



% = 



l + i\u2\ iu*2 
-iu2 l-i\u2\ 



Ml = \U2\, 



(A6) 



cosh \u2\ 
ie"^sinh|M2| cosh 



sinh \ u2 



U2\ 



ui = 



with e*'^ = U2/\U2\. 

Let us also consider two further limiting cases, when 
the range of the forward scattering potential ui{x) is 
much longer (shorter) than the range of the backward 
scattering potential ^2(2^), i-c ai ^ 02 (ai <C 02 ). In 
these two cases we get 



Tn = 



e*"icosh|u2| ie "^sinh|it2| 
-ie*"^ sinh \u2 1 e"™^ cosh \u2 



(AT) 



for fli ^ 02 and 

J- _ ^ /^cosh |u2| -I- i tanui ie~*"* sinh |u2| 
" "'^ \^ — ie*'^ sinh |u2| cosh |u2| — i tanui 

(A8) 

for ai <ti a2- By construction, the determinant of the 
transfer matrix is unity for arbitrary Mi,2(a^)- 

The elements of the scattering matrix for the incident 



waves 



^ and 



at a: < and a; > 0, respec- 



tively, are defined by the equation 



(A9) 
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Hence the scattering matrix is determined by the ele- 
ments of the transfer-matrix T as 



S: 



t f 
r i 
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1 Ti2 
-T21 1 



In the above two cases with ai = 02 we have 
1 I I iu 



S = 
S = 



1 - i\u2\ V*"2 1 

l/cosh|it2| ie"**^ tanh |u2 
ie*"^ tanh |u2| l/cosh|u2| 



ui = |U2|, (AlO) 

Ml = 0. 



Whereas for model potentials of different range we obtain 



S 



1/ cosh |lt2| 
ie'-^tanhlual 



ie *'^tanh|M2 
1/ cosh \U2\ 



(All) 



for ai ^ a2 and 



S 



1 



1/ cos Ml 
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sinh |m2 



cosh |m2| - ^tanui ^ie'* sinh |u2| 1/cosmi 

(A12) 

for ai <C a2- We see that if the forward scattering mi 
has wider range, than M2, then Eq. (|All|) differs from 
the second Eq. (jAlOp only by the overall phase which is 
unimportant in physical observables. At the same time, 
if the model regularization of ui is narrower than U2 , the 
resulting Eq. (jA12p closely resembles the S'-matrix in the 
double-barrier situationJ^ 

Expanding the reflection coefhcient in powers of ui_2 
and using ((X5|l . (|MT|) . (^12)) we have 



\U2\ - g|U2| 



ai » a2. 



|r| ~ \u2\- -\u2\iul + 2\u2\^), ai = a2, (A13) 
6 

|r| ~ \u2\-l\u2\{3ul + 2\u2\^), ai<a2, 
6 

which shows that the forward scattering amplitude ui 
enters the reflection coefficient in the third order, i.e. 
beyond the Born approximation. Moreover the contri- 
bution of Ml to the S'-matrix depends also on the details 
of the microscopic Hamiltonian and the regularization of 
the model (5-function potentials. This complication hap- 
pens at the step (IA2[) : however, once the transfer matrix 
Tq in (jA4[) is known, the connection to the scattering 
state basis is straightforward. 

Returning to (g]), we notice that "01(2;) obeys 

the equation of motion for right movers, and the spinor 
^ft^ = (■0+(a;),r-0+(-x)) is an eigenstate ((X2|) of (|XT|) 
with the above property (|A9p . We thus clarified the scat- 
tering states representation in the framework of the chiral 
Hamiltonian ([T]). 

In the bosonization formulation, ([5]), the assumption 
is made that the forward scattering Mi(a;) can be com- 
pletely removed from the Hamiltonian. This statement 
is not obvious, in general. Consider the redefinition 
f{x) = tf{x) — uisgnai{x)/2, which eliminates the term 



with Ml from Eq. ([5]), where the regularized sign function 
sguaiix) = — 1 -|- 2 J^^dy Sai(y) is used. We have the 
backward scattering in the form 

M2(a;) cos(2iy9(a;)) = M2(a;) cos(2(^(x) — MiS(7nai(a;)), 
U2 5{x) cos{2(p{x)), 
— > M2 (5(a;) cos(2(^(x)) cos Ml, (A14) 

for the above cases of regularization with ai ^ 02 and 
ai ^ 02, respectively. 

APPENDIX B: UNITARY TRANSFORMATION 

We consider a transformation of the operators Jj 
Oj{x) = UUj{x)U 

with 



U = expi / dx dy 0{x — y)Bj{y)Jj{x), (Bl) 

J — oc 

/oc 
dyB,{y)e,{y), (B2) 
-00 

$,(2/) = / dxJ^ix). (B3) 

In our case we have Bj{x) = Bj5{x) and we may denote 
the direction of B as "1", i.e. B||ei. Then we have the 
following set of equalities from the Kac-Moody relations 

[Ji(y),$i(a:)] = ^5{x - y) 
[Ji(y),z j dxB,{x)Q,{x)] = -j^B,{y) (B4) 

My)U - U{My)-^B,{y)) 

And for J^{x) = J3{x) + iJ2{x) one has 

[J+{x),My)] = J+{x)S{x~y) 
[J+{x),ei{y)] = J+{x)d{x-y) (B5) 

[J+{x),i I dy Bi{y)ei{y)] = tj+{x) f dy B^{y) 



So that the relations between the operators are 

OAx) = UUi(x)U = Ji(x) - —BUx) (B6) 

47r 

03{x)+i02ix) = Ulj+{x)U (B7) 
= {J3{x) +iJ2{x))eyipi dy Bi{y) 



The transformation U shifts the value of Oi{x) by a 
c— number, so that strictly speaking, the operators Oj(x) 
do not obey the same Kac-Moody algebra in the region 
where Bi{x) 7^ 0. It is clear, though that this might af- 
fect only the physically non-observable component of the 
fermionic density and is therefore not important. 
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APPENDIX C: EQUATIONS OF MOTION IN 
CURRENT OPERATOR FORMALISM 



In this section we Jbcus on the time evolution of the 
observables J3{—x), J3{x), Eqs. P^ . P5|) . using the rela- 
tion idtA = [A, H]. It is more convenient to use the rep- 
resentation dlZl), with 0{x) = J3(x)e(-a;) + J3(a;)e(a;). 
In terms of operators Oj (x) one has 



Ps{x) = (Osix) - Osi-x)) 
js{x) = vf{03{x)+03{-x)) 



(CI) 



Using (I20|) and formulas of Appendix it can be easily 
shown that we have —dtOslx) — dx03{x) in the leads, 
|a;| > L, whereas in the interacting region a < |a;| < i 
we obtain 



dtOsix) = dMx) ~ gdM~x), 



(C2) 



where we let vp = 1- Since 03(2;) is the observable den- 
sity, the relation (|C2[) looks deceptively close to the de- 
sired answer. However, in the region of the barrier |x| < a 
(with B||Oi) we have 

^dtO^ix) = dM^) + B{x)o^(x), 

and the O2 component comes into play. (In other words, 
for the partially reflecting barrier, the observable of our 
problem, Jz(x)^ includes the eigenmode operator J2{x).) 
Therefore it is necessary to consider also the time evolu- 
tion of operator O2, which reads 

- 9*02(2;) = 9x02 (x) - 4^.gOi(a;)03(-x), (C3) 

at a < < L and —dt02{x) ~ dx02{x)A-B{x)0z{x), at 
|a;| < a. The four-fermion object Oi(x)03(— x) is not re- 
duced to 02{—x) and therefore we do not obtain a closure 
of the equations of motion. 

Our situation is somewhat similar to the well-known 
Heisenberg model of spins with pairwise exchange inter- 
action. Due to the non-commutative character of spin 
operators, the attempt to write the exact time evolution 
of spin leads to an infinite series of coupled equations, 
including ever increasing numbers of spins. The approx- 
imate solutions are obtained by truncating this series at 
some step, according to a certain criterion, and discard- 
ing higher terms as inessential)^ In magnetism, these 
criteria might be (i) low temperatures in the magneti- 
cally ordered phase, T <C Tc, (ii) large value of spin S or 
(iii) wide range of interaction. In our case, we choose the 
criterion of relevancy of operators, i.e. their importance 
in the logarithmically divergent theory. 

We now briefly discuss the dynamics of 04(2;) = 
47rOi(x)03(— x). The compound operator 04(x) con- 
sists of four fermion operators and may be decomposed 
into the normally ordered part and a part, representing 
internal contractions. To clarify the notation we write 
■01(2;) = a|(x), ip2{—x) — ai{x), then we have 

04(2;) — 7r(a|(.T)a|(x) + a|(2;)a|(2;)) 

•K{a\{~x)a^{-x) - a\{-x)ai{-x)) (C4) 



Contracting here, we get 04(x)— : 04(2;):= —x ^02{x) 
with (cf. Eq. (EH)) 

02(x) = ^\a\{x)ai{-x) + a\{-x)ai{x)) + h.c. (C5) 

which does not coincide with 02(2;) ~ {2i)^^a}-^(x)a\^(x)-\- 
h.c. We may write 



02(X) ~ 02(-x) - xdM-x) + 



(C6) 



and, dropping the less relevant higher derivatives, obtain 
from CI 



- 9*02(2:) ~ dMx) + g{x-^ - dx)02{-x) ~g:Oi{x): 

which means that the most relevant part of the operator 
04(2;) has the form 



04(2:) ~ {~l/x + dx)02{-x). 



(C7) 



Regarding the normally ordered part : O4 (x) : , its time 
evolution is given by 

-dt-.Oiix): = 4^:03(-2;)9,Oi(2;)-Oi(2;)9,03(-x): 
+ATrg 105(2;): 
05(x) = A7:Oiix)Ol{-x) + Oi{x)dMx) 

i.e. described by new operators, not reduced to simple 
spatial derivatives of : 04(x) :. Such derivative objects (i) 
obey an infinite set of coupled equations and (ii) possess 
formal scaling dimension greater than two, and hence are 
irrelevant in the RG sense. 

Our strategy of such truncation (|C7p is successful only 
in part, as we explain in the next Appendix [P] 



APPENDIX D: SOLUTION OF THE EQUATION 
OF MOTION 



Taking a Fourier transform in time t, and introducing 
the vector 

^'(X,w) = [03{x,Uj),03{-X,U!),02{x,U!),02i-X,Uj)] 

the above truncated set of equations can be written as 

dx'^ix^Lu) = Do{uj,x)\g=o'i'{x,Lu), \x\>L (Dl) 
i)i9:r*(x, uj) = Doi^^, x)^{x, Lu), fl < |2;| < L 

with 



Do{ll>,x) 



iLu, 0, 0, 

0, ~iuj, 0, 

0, 0, iuj, ~gf{x) 

0, 0, ~gf{x), -iLu 

1, -g, 0, 
-g, 1, 0, 

0, 0, 1, ~g 



1 u, 

V 0, 



(D2) 



(D3) 



0, -5, 1 
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and the fermionic correlation function at finite tempera- 
ture 



2t:T 



(D4) 



sinh 2nTx 

The solution to (|D1[) is given by the x— ordered exponent 

^ix,u;) = U{x,y)^iy,u) (D5) 
Uix,y) = T.exp [ dz D^^ bo{uj, z), (D6) 



The overall transfer-matrix over the interaction region is 

U{L, -i) = ij{L, a).U{a, -a).U{^a, ~L) (D8) 

Finally we observe that the definition of ^'(x) implies that 
the transfer matrix connecting opposite points should 
satisfy the boundary condition: 



U{x,-x)^{-x,uo) = /.*(-a;,w) 



similarly to (|A2[) . but now with transfer-matrix U{x,y) 
for densities. The transfer-matrix in the region of the 
barrier < a is already known and given by 



/O, 1, 0, ON 

1, 0, 0, 

0, 0, 0, 1 

Vo, 0, 1, oj 



(D9) 



(DIG) 



U{a, —a) 



/ cos 26*, 0, sin 26*, 

0, cos 26*, 0, -sin 261 

-sin 26*, 0, cos 261, 

V 0, sin 26*, 0, cos 26* 



(D7) 



and we should choose a solution, by demanding that the 
scattered component, O2, is absent in the incoming wave, 
i.e. ^'3(— cx), Lu) = 0. 



J 



We consider here only the technically feasible limit, uj 
X > L. In the interacting region 



U{L,a) 



/I, 0, 

0, 1, 



0, 
0, 



0. In this case U{x, L) 
\ 



U{—L, —x) = 1 in the leads. 



0, 0, e^f cosh gA, -e'^* sinhgA 
Vo, 0, -e-s'^sinhc/A, e-^'^coshgA / 

^ , „, , , 9s 1 , / tanhTrTiN 



(Dll) 



(D12) 



A similar contribution on the negative semiaxis, U(—a, — £), is obtained from (jPlip by the change A —A. Overall 
we have 



U{oo, —00) = 



( cos 26*, 

0, 

e^f ^ sin 20 cosh 5A, 



0, 
cos 20, 
e~9^^sin26'sinh5A, 



sm 261 cosh 5A, -ea'-i^ sin 26* sinh gh\ 
sin 26* sinh 5A, e^^^^ sin 26* cosh 5 A 



-e a "sin 20 sinh gA, — e gj^^ 20 cosh gA, 



V-e"f ^sin20sinhgA, -e^^ A 
Solving Eq. (jD9p together with the above condition 5'3(— cx)) = 0, we find 

5-(-oo) = (cos^ -t- e^^^ sin^ 0, cos^ - 
which is the desired solution. Let us briefly discuss it. ^ 



cos 20, 
0, 





cos 20 



2s^sin2 0,O, 



.9(1-5)A 



sin 20) 



(D13) 



We observe that the d.c. current ji."^) — 0-i(x) + 
0-i{—x) is independent of a;, as it should be. This prop- 
erty follows from the form of matrix (|D11|) and condition 
(jPgp . The relation of the current to the voltage bias de- 
fines the conductance of the system. In our definitions, 
the voltage bias is proportional to the difference in num- 
ber of incoming left- and right-going electrons, which is 



203(— 00) — 2'I'i(— cxd), so that the conductance is 

^ 3{x) ^ 1 / 03(^3) \ 
203(-oo) 2 \0^{-oo) ^ ) 



cos2 e29A sin^ ^ ' 

which coincides with the above result (fTTj) . 

At first glance, the result (jD13p should thus be viewed 
as satisfactory. However there are details, involved in 
its derivation, which spoil the significance of this ap- 
proach. First we notice that, by keeping less relevant 
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gradient term in the above relation, (jC7p . we get the ap- 
pearance of terms in the definition of A, eq. (|D12p . 

and also in the exponents, e^^^, in (jDll[) - (|D13[) . Fur- 
ther, trying a different approximation instead of (jC7[) . 
e.g. 04(3::) — {—l/x + dx)02{x), would not produce the 
above renormalized value of conductance. 

Therefore we conclude that using the equations of mo- 
tion method and persisting with the current operator al- 
gebra leads to ambiguities, connected with the non-local 
character of the interaction. 



APPENDIX E: CALCULATION OF 
CORRECTIONS: ZERO TEMPERATURE 



Proceeding as described in Sec. IV Al we find correc- 
tions to the conductance, or to Y in (PT|) . In the first 
order of g the correction from the diagrams of the type 
{4} has the form 



dx 

X 



^g(l - r2)Ao (El) 



processor), we obtain 



1 



{xi+X2+X-i,Y {Xl+X2-X^Y 

4 



-g3(i-r^)(i + r2)A£, 



(E4) 



-g\l-Y^f jWdx, 



1 



Xi{x2 + X:i){xi + X2 + Xs) 



(symm.) 



,1 



-g3(i_y2)2_(A2 + 2Ao(l-ln2)), (E5) 



, (symm 



1 



Xi{x2 + X3Y 

= g\l-Y^)Y^Ka{^a-2\n2), (E6) 



In the second order of g we have two sets of diagrams 
for corrections to the conductance. The diagrams con- 
taining two fermionic loops lead to the following expres- 
sion : 



^ g^Y{l-Y^) I 



dxidx2 

+ X2f 



-g^Y{l - Y^){Ao - 2 In 2) 



(E2) 



where xi^2 are two points of fermionic interaction. 

All diagrams consisting of only one loop, {6} in above 
notation, produce the correction 



SYi^^ = -c 



'-Y{l-Y^) f 

J a 



dxidx2 
X1X2 



-g'Y{l - Y^)Al 
(E3) 

We should mention that individual diagrams in the set 
{6} (and higher orders) may contain singularities of the 
form (xi — 2^2)^^ etc. However these singularities cancel 
each other in the resulting expressions. For example, one 
finds two particular contributions Xi^{xi — X2)~^ and 

X2^{X2 — Xi) 

form, {X1X2) m 

In the third order, we have to analyze the following 
types of the diagrams : {4,2,2}, {4,4}, {6,2}, {8}. We 
assume the ordered sequence, L > xi > X2 > X3 > a, 
in all integrals below. Performing rather long computer 
calculations (the stable routine in Mathematica requires 
about ten hours of computation time for dual-core 3 GHz 



whose sum leads to the above "regular" 
^-1 (El. 



<5rW = 2g\l-Y^) j \{dx. 



i-3y2 



XlX2X'i 



l-Y' 



{Xi + X2){X2 + X'i)(X'i + Xl) 

1-3^2 



(E7) 



g^l-Y') 



In the final expressions (jE4[) . (jE5[) . (jE7[) we omitted finite 
parts, 0{K^), which are unimportant in the order g'^. 
These expressions lead to Eq. (^5]) above. 



APPENDIX F: CALCULATION OF 
CORRECTIONS: FINITE TEMPERATURE 

In case of finite temperature we proceed in the same 
way as for T = , only with the integration replaced by a 
summation over Matsubara frequencies. Introducing the 
scaling variable A according to (|47| we have for a ~ 



^^{4,2} ^ g^Y{l-Y^) I 

Ja/^ sinh 2(.xi -t- X2) 

= ]^g^Y{l-Y^)(K-i\n2), 



(Fl) 



where Xi^2 are two points of fermionic interaction, and 
Xi = Xiji^. In the last equality in the right-hand side of 
Eq. (jFl[) . we shall neglect the terms ~ 0(1). 
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Instead of Aq in the expression (p3l) for SY^^^ we now 
have the foUowing integral 



cosh(a;i — X2) 



^'^ ,_ ,_ 3cosh(xi +^2) 
jj/^ sinh(2.Ti) sinh(2.f2) cosh(.Ti + X2) 

dTldT2 { ^ 1 



2 



1 



T1T2 

2 1 + T1T2 



(F2) 



= A^-F[Ti,ro] 



with Ti — tanh^i . The integration hmits in (|F2|) 
are tq = tanh(a/^) and = tanh(L/^). Here 
F\tlM =^L2{-totl) - (i2(-T|) + L2(-To2))/2 and 
L2{x) = is the dilogarithm function. We have 

F[tl,to] = in the hmit T < vp/L. 

In contrast to the T = case, the expressions now be- 
come extremely cumbersome in the third order of g. The 
summation over Matsubara frequencies lead to expres- 
sions involving sums and products of exponents exp(iEi). 
A general criterion for the correctness of intermediate 
expressions is the property a change of sign of each indi- 
vidual diagram upon performing the reflection transfor- 
mation Zi — > —Zi. 

An additional check of the validity of calculation is the 
observation that we should get zero correction to G = 
1 in the absence of the barrier, 9 — 0, and to G = 
for the fully reflecting barrier, 9 = 7r/2, for arbitrary T. 
This is not a trivial statement, because each diagram in 
third order contains a prefactor, either cos 86', or cos46' 
or cos 00 = 1. We verified that the algebraic expressions 
canceled each other in this case. 

The actual expressions are not shown here, and we 
describe our method of analysis below. 

Let us consider an expression of the form of a triple 
integral, relevant to our discussion in Sec. |V] 



Therefore the recipe for the determination of c is first 
to determine the limiting form of /2, taking Z2 ~ a ^ ^ 
; in this case the expressions are drastically simplified. 
This will determine c at T — > 0. To this term we should 
add the contribution lima^o o,f[a, Z2, 23], integrated over 
Z2 in the interval (0, 00) ; the sum of these two contribu- 
tions gives the value of c at finite temperatures, L 3> 



APPENDIX G: "NON-UNIVERSALITY OF 
HIGHER TERMS" IN /3-FUNCTION 

In a very general sense the /3-function of the RG- 
equation is non-universal: it depends on the precise def- 
inition of the scaling quantity g . Let the RG equation 
read 

dgr/dA = (3{gr) = 62.9^ + hdr + hgt + ■■■ (Gl) 

and consider a change of variable gr — 9 + C2<f' + c^g^. 
With the same accuracy we have 

dg/dA = Pig) = P{g)/{d9/dg) - 623' + hf 

+ ih + hC2+b2{cl-C3))f + ... (G2) 

which shows that only the two first two coefficients &2 and 
63 are invariant. Since 64 is associated with the three- loop 
contribution, one may speak about a non-universality of 
the /3-function beyond two-loop order. 

Let us elucidate the origin of the Eq. (jG2[) . and 
compare it with the equation obtained in the Callan- 
Symanzik approach. In the Gell-Mann-Low approach 
we perform our calculations with the running coupling 
constant gr and add counterterms to the Hamiltonian 
in order to achieve the finite value of the bare constant 
go- One can show that the above value of dgr/dA = 
~{dgo/dA)/{dgo/dgr) is obtained from the expression 



/(«) 



dzi 



dz2 / dz3f[zi,Z2,Z3]. (F3) 

Zi J Z2 



We assume that /(a) ~ — clna + . . . at a — s- and wish to 
determine the coefficient c in this asymptotic behavior. 
We have 



df{a) 



a^o din a 



lim 



lim 



dZ2 



dzs af[a, Z2, Z3] 



. (F4) 

Here taking the limit in the integrand may lead to incor- 
rect results. We have two possibilities for the behavior of 
the quantity /2[a,Z2] = J^^^z^afla, Z2, z^], for small a, 
namely f2[a 0, Z2] ~ a/{z2 + a)^ and f2[a 0, Z2] ~ 1. 
The first type of contribution, /2 ~ a/{z2 + a)^, is al- 
ways important and survives in the limit of T 0. 
Recalling that /[a, 22,2:3] is exponentially suppressed at 
^ see that the case f2[a,Z2] ~ 1 contributes to 

c only at finite temperatures, i 3> ^ ; typically, we find 
f2[a^0,Z2] ^ (ecosh(2z2/C))-^ 



go = 9r - h2Agl - {bsA - blA^)gl 



-ibiA--b2b3A^ 



hlA^)gt 



(G3) 



If we make the above change g^ ^ g the right-hand 
side of (jG3ll and use the same recipe for dg/dA, we return 
to (lG2l) . 

In the CS approach we fix the bare value go and calcu- 
late the renormalized quantity gr{go,A). Inverting the 
equation gQ{gr,A) and using the same recipe we find 
dgr/dA. One can check that 



9r = go + b2Agl + {b3A + blA'^)g^ 
+ {b4A+h2b3A^ + blA')g^- 



(G4) 



Assume now that we have changed the starting value 
9o — 9Q+C29o~^C3gQ + . . ., according to above prescription. 
Evidently, the expansion for gr{go,A) will be different 
from (jG4p . as well as the inverse function go(ffr,A) will 
not coincide with (|G3|) . One can verify however that the 
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/?-function, —{dgQ/dA)/{dgQ/dgr), is given by the same 
Eq. (|G1[) and not by Eq. (|G2[) . This fact becomes rather 
obvious, when we notice that only the left-hand side of 
(|G3[) is changed and it corresponds to a change in the 
boundary condition for (7r(A) rather than to a change in 
the differential equation for it. 



APPENDIX H: COMPARISON TO CFT 
SOLUTION 

Consider the differential equation 

[-d^ + se^y + ey]^{y)^0 (HI) 

such that ^'(y oo) = and \E'(?/ — oo) — y ~ yo- 
Here 

s = X2exp(2(l - K)A) 

and A = ln(To/T). For given s one finds a solution ^'(y) 
characterized by a unique value yo- Then, according to 
Lukyanov,- - the conductance can be determined as 



G{T) = 1 + Ks 



dyo 
ds 



1 



K 



dyo 



2(1 -X) dA 



(H2) 



We now fix the overall temperature scale r = T/Tq 
by the high-T behavior, G ~ 1 — t^(^~^\ and ask what 
is the coefficient C in the asymptotic low temperature 
expression G ~ Ct"^^^ To answer this question, we 

slightly rewrite the differential equation. Introducing a 
new variable z = 2e^/^ we have 

[-^2 - z-^d, + 1 + s(z/2)2(^-i)]*(z) = (H3) 

Ats = (i.e. T — > oo) the solution is a modified 
Bessel function ^(z) = —2Ko{z). At infinitesimal s 



we seek a solution in perturbation theory as ^{z) = 
—2Ko{z) + sR{z) + 0{s'^). After some calculation we find 
the conductance G(T) ~ l-KsJ^{K) with J^{K) = i?(0) 
given by 

f°° riK 

HK) - T^[K] / ^F^{K,K,l--n)-— 
Jo 1 + 

= T^[K]/T[2K] . (H4) 

In the opposite limit, s ^ oo (T = 0), we choose a 
variable z = {2^/s/K)e^y/^ and arrive at an equation 



+ l + si(z/2)2(^"'-i)]^(z) = 0, (H5) 



with si = s-i/^7;i'2(if-i-i) ^ 2. The solution is found 
similarly to ([HSj and we get G{T) ~ -sdsSiJ^iR-^). 
Comparing the two solutions we obtain the prefactor 



(H6) 



In the particular case of ii" = 1/2 we have C = 7r'*/48 = 
2.029..., in accordance with Ref. For K = 1/3 

we get C — 10.064558 . . ., a value slightly different from 
C = 10.0638 reported inM^, where it was obtained by 
numerical solution of an integral equation. If we take 
cg = 1/4 in dn?]), then we would obtain G = 2 (G = 9) for 
the cases of iiT = 1/2 {K = 1/3), in good agreement with 
the above CFT values. It should be noted, however, that 
apart from this agreement in the matching coefficient G, 
the theory in^ provides the value of unity for G in the 
clean wire, whereas the previous work o^^i^^'^^ gave G = 
K in that case. 
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